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A theory of crossed Andreev reflection in structures consisting of a superconductor with two 
ferromagnetic leads is presented. The electric current due to the crossed Andreev reflection strongly 
depends on the relative orientation of the magnetization of two ferromagnetic leads. It is shown 
that the dependence of the electric current and the magnetoresistance on the distance between two 
ferromagnetic leads is understood by considering the interference between the wave functions in 
ferromagnets. The current and the magnetoresistance are calculated as functions of the exchange 
field and the height of the interfacial barriers. 
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I. INTRODUCTION 

Much attention has been focused on the spin depen- 
dent transport through magnetic nanostructuresii The 
tunnel magnetoresistance (TMR) was observed in fer- 
romagnet/ferromagnet (FM/FM) tunnel junctionsSiSi^. 
In ferromagnet/superconductor (FM/SC) tunnel .junc- 
tions, the current flowing thorough the tunnel junctions is 
spin polarizedf£ When the spin polarized quasiparticles 
(QPs) is injected into SC from FM, the superconduct- 
ing gap is suppressed due to the spin accumulation in 
FM/SC and FM/SC/FM j unctions -ii2i2iiaJii2*iii4 The 
detail studies of the spin transport and relaxation in SC 
have been doneii^iiSiii 

In recent years, many theoretical and experimental 
studies in relation to Andreev reflectioniS in FM/SC 
metallic contacts have been done because the spin polar- 
ization of conduction electrons is estimated by measuring 
the conductance in this system,i9i2°i'^b22i23,24,25,26,27,28,29 

In FM/SC/FM double junction systems, the coherence 
length in SC is extracted by measuring the magnetore- 
sistance (MR)ii22i^ In a system consisting of SC with two 
ferromagnetic leads FMl and FM2 (see Fig. there is a 
novel quantum phenomenon called the crossed Andreev 
reflection|32,33,34,35,36,37,38,39,4Q,4i,42,43 ^^^ven an electron 

with energy below the superconducting gap in FMl is in- 
jected into SC, the electron captures an electron in FM2 
to form a Cooper pair in SC. As a result, a hole is cre- 
ated in FM2. Deutscher and Feinberg''- have discussed 
the crossed Andreev reflection and MR by using the the- 
ory by Blonder, Tinkham, and Klapwijk (BTK).^^ They 
argued that the crossed Andreev reflection should oc- 
cur when the distance between FMl and FM2 is of the 
order of or less than the size of the Cooper pairs (the 
coherence length), and calculated the probability of the 
crossed Andreev reflection in the case that both ferro- 
magnetic leads are half metals and the spatial separation 
of FMl and FM2 is neglected (one dimensional model), 
i.e., the effect of the distance between two ferromagnetic 
leads on the crossed Andreev reflection is not incorpo- 
rated. Subsequently, Falci et alfi^ have discussed the 



crossed Andreev reflection and the elastic cotunneling in 
the tunneling limit by using the lowest order perturbation 
of the tunneling Hamiltonian. However, to elucidate the 
effect of the crossed Andreev reflection on the spin trans- 
port more precisely, it is important to explore how the 
crossed Andreev reflection depends on the distance be- 
tween two ferromagnetic leads as well as on the exchange 
field of FMl and FM2, for arbitrary transparency of the 
interface from the metallic limit to the tunneling limit. 

In the present paper, we present a theory of the crossed 
Andreev reflection in structures consisting of SC with two 
ferromagnetic leads. By extending the BTK theory to 
this system, we derive an expression of the electric cur- 
rent and calculate the current and MR originated from 
the crossed Andreev reflection. The dependence of the 
current and MR on the distance (L) between FMl and 
FM2 is examined. It is shown that the dependence of the 
crossed Andreev reflection on the distance L comes from 
the interference between the wave functions in FMl and 
FM2, and the probability decreases rapidly as [kpL)'^ 
with increasing kpL, but not the coherence length of 
SCpS, where kp is the Fermi wave number. The current 
and MR are calculated as functions of the exchange field 
and the height of the interfacial barriers in order to clar- 
ify the crossed Andreev reflection in the spin transport 
of the present system. 



II. MODEL AND FORMULATION 

We consider a system consisting of a superconductor 
(SC) with two ferromagnetic leads (FMl and FM2) as 
shown in Fig. [T] FMl and FM2 with width Wp are 
connected to SC with width Ws at a: = 0. The distance 
between FMl and FM2 is L. The system we consider is 
described by the following Bogoliubov-de Gennes (BdG) 
equationii^ 
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FIG. 1: Schematic diagram of a superconductor (SC) with 
two ferromagnetic leads (FMl and FM2). FMl and FM2 
with width Wf are connected to SC with width Ws at a; = 0. 
The distance between FMl and FM2 is L. 



where Hq = — (?i^/2m)V^ — is the single particle 
Hamiltonian measured from the Fermi energy fip, E is 
the QP excitation energy, and a = +{—) is for the up 
(down) spin band. The exchange field /lea; is given by 



(a;<0, \y~L/2\<WF/2), 
(x>0), (2) 
{x<0, \y + L/2\<WF/2), 



where +ho and —ho represent the exchange fields in FM2 
for the parallel and antiparallel alignments of the mag- 
netizations, respectively. The superconducting gap is ex- 
pressed as 




A(r) 



(x>0, |y| < W^5/2) 
(x < 0) . 



(3) 



We assume that the temperature dependence 
of the superconducting gap is given by A = 

Aotanh (l.7A^/Tc/T - where Aq is the supercon- 

ducting gap at T = and Tc is the superconducting 
critical temperature. In order to capture the effect of the 
interfacial scattering, we employ the following potential 
at the interfaces, a; = 0: 



HB{r)^H6{x){9i{y) + 02{y)}, 



(4) 



where d{x) is the delta function and di(2){y) = 
9{WF/2-\y- {+)L/2\), 6{x) being the step func- 
tion. Throughout this paper, we neglect the impu- 
rity scattering in SC and the proximity effect near the 
interfaces,33i47,48,49,50,5i 

The solution of the BdG equation in the SC region is 
given by 



*±,+ (r)= ("")e±^'=^<I>sc,Ky), 



(5) 



where r = [x, y), and uq and uq are the coherence factors, 

(6) 
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For i? < A, Mo and vq are complex conjugates. $sc,;(2/) 
is the wave function in the y direction. 
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where I is the quantum number which defines the channel. 
The eigenvalue of the y mode for channel I is 



E, - 



2m \Ws 



(8) 



The x component of the wave number of an electron 
(hole) like QP, k^''- \ is expressed as 

/2m 



±\/E^ - - El. 



(9) 



In the FMl (FM2) region, the solutions are given by 

= ±'<i^ $FM1(FM2), 
3±^P^,i^$pMi(FM2),i(2/), 



(r) = 



(10) 



where $ 
tion 



FMl(FM2),i 



(y) is the wave function in the y direc- 



$FMl(FM2),i(y) = ^/-J^sin-i^ 



and ' is the x component of the wave number of an 
electron (hole) with a spin; 
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We consider the scattering of an electron with a spin 
in channel n injected into SC from FMl. There are the 
following six processes: the ordinary Andreev reflection 
and the normal reflection at the interface of FMl/SC, the 
crossed Andreev reflection, the crossed normal reflection, 
the transmission to SC as an electron like QP, and the 
one as a hole like QP. Therefore, the wave function in 
each region is expressed as follows: In the FMl region. 



*FMi(r) 



J j e<"- FMl Ay) 
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+ b^,i, 



$FMl,i(y), (13) 
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in the FM2 region, 

oo 

^'FM2(r) = ^ 



and in the SC region, 

oo 

«'sc(r) = 

1=1 L 



$FM2,i(y), (14) 



^scAv)- (15) 



Let us evahiate the current in FMl. When the bias 
voltage V is apphed to the system (see Fig. 0, the cur- 
rent carried by electrons with a spin in channel m is given 
by 



l,e 



{E)]dE, (19) 



where h is Planck constant, and /^'^j _^ (E) is the distri- 
bution function of an electron with positive group veloc- 
ity in the x direction and is expressed as 



l,e 



(E) = h{E- eV) 



(20) 



where /o (E) is the Fermi distribution function. The dis- 
tribution function of electron with negative group veloc- 
ity in the x direction /j:',^ ^ (E) is given by 



Here, p^^, q^^, and are the wave numbers in FMl, 
FM2, anii SC', respectively. 

The boundary conditions at the interfaces {x — 0) are 
as follows: 



*FMi6'i(y) + *FM26'2(2/) = "^scOsiy), 

'^Siy) - -r [*FMiei(2/) + ^FM2^2(2;)] 



(16) 



dx 



dx 



[*FMl6ll(y)-|-*FM26'2(y)], (17) 



where Os{y) ~ {Ws/2 — \y\). From the boundary con- 
ditions, the coefficients aa.in, ba,in, Ca,in, dcr,in, Qio-.in, and 
Pa,in are determined (see Appendix)^24^2iS4^. The prob- 
abilities of the Andreev reflection R]^''^^, the normal re- 
flection i?;^',mnj the crossed Andreev reflection i?;^;TO„, the 
crossed normal reflection R];^^^, the transmission to SC 
as an electron like QP, Tj^'mn' ^'^'^ t;he one as a hole like 
QP, T'jJmnj are written as. 
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where vs^i and Vp ; are the group velocity of an electron 
in channel / in SC and the one with a spin in channel 
I in FMl, respectively, Ns,i and Npi are the density of 
states in channel I in SC and the one of a spin band in 
channel I in FMl, respectively. Using the relations. 
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and the conservation law of the probability. 
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where the superscript e'{h') and 1 in Eq. (|18|l indicate 
the electron (hole) like QP in SC and the injection from 
FMl, respectively. 
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we obtain 
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The current carried by holes with a spin in channel m in 
FMl, /^'^j, the currents carried by electrons and holes in 
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FM2, /^f„ and I^'l^^, respectively, are calculated in the 
similar way as 
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By using Eqs. (|24|I - H27|I . we obtain the total current in 
the system 



j^y- rjl.e _^j2,e _^j2,hl 

We define the magnetoresistance (MR) as 
^AP — Rp Ip — Iap 



MR; 



Rp 



(28) 



(29) 



where i?p(AP) — ^/^p(ap) is the resistance in the parallel 
(antiparallel) alignment. 



III. RESULTS 

In the following calculation, we take the tempera- 
ture, the applied bias voltage, the width of SC, and the 
superconducting order parameter to be T/Tc = 0.01, 
eV/Ao = 0.01, Ws = 1000/fcF, and /if/Ao = 200, re- 
spectively, where kp is the Fermi wave number. First, 
we consider the case that FMl and FM2 are half metals 
{ho/^ip = 1) a-nd the strength of the interfacial barrier 
Z = mH/h^kp = 0. The width of FMl and FM2 is 
taken to be Wp ~ A/kp. In this case, there is only one 
propagating mode (^ = 1 in Eq. (|10|l ). We obtain the 
maximum possible value of MR, i.e., MR = — 1 indepen- 
dently of L. In order to understand this behavior, we 
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FIG. 2: The current as a function of L. FMl and FM2 are 
half metals {ho/fJ-p ~ 1). The solid and dashed lines are for 
the currents in the antiparallel and parallel alignments of the 
magnetizations, respectively. 



consider the L dependence of the currents in the paral- 
lel and antiparallel alignments as shown in Fig|21 When 
an electron with up spin in FMl is injected into SC, the 
ordinary Andreev reflection does not occur because elec- 
trons with down spin are absent in FMl. In the parallel 
alignment, the crossed Andreev reflection does not occur 
either because there are no electrons with down spin in 
FM2. Therefore, no current flows in the system as shown 
in Fig 121 On the other hand, in the antiparallel align- 
ment, while the ordinary Andreev reflection is absent, 
the crossed Andreev reflection occurs because there are 
electrons with down spin in FM2, which is a member of 
a Cooper pair, for an incident electron with up spin from 
FMl, and therefore finite current fiows in the system as 
shown in Fig [21 As a result, we find MR = — 1 irrespec- 
tive of L in the case of half metallic FMl and FM2. The 
current in the antiparallel alignment decreases oscillating 
with increasing L. The behavior of the current is under- 
stood as follows. From Eqs. H24|) - (|28|l . the current in 
the antiparallel alignment at low temperatures and low 
applied bias voltage is expressed as 



/af 



62, /le 
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2,eh 
T,ll 



(30) 



It is shown from Eq. (|A1|) - (|A16|I in the Appendix that 
the L dependence of the probability of the crossed An- 
dreev reflection R^^'^l is originated from the interference 
term between the wave functions of FMl and FM2, and 



R^'^li (X [1 



sm{2kpL + (f))]{kpL) cxp(-L/0, (31) 

where ^ = ^GL(7rA/2%/A2 - E'^), ^gl = hvp/i:A being 
the Ginzburg-Landau (GL) coherence length, vp is the 
Fermi velocity, and is a phase defined as Eq. ljA16|l . 
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FIG. 3: The absolute value of MR as a function of L in the 
case that the exchange field ho/i-iF are 0.4, 0.5, 0.6, 0.7, 0.8, 
and 0.9. 
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FIG. 4: The current as a function of L in the case oiho/fiF ~ 
0.6. The solid and dashed lines are for the currents in the 
antiparallel and parallel alignments, respectively. 

L dependence as R^'ii in Eq. (|81(l . and therefore the cur- 
rent in the antiparallel alignment H3U|I decreases rapidly 
with a rate of {kFL)~^ oscillating with period of tt with 
increasing kpL- Note that the kpL dependence of the 
probabilities is dominated by the term {kpL)'^ , not the 
exponential term cxp(— L/^),-^ since kpS, ^ 

We next consider the L dependence of MR for several 
values of the exchange field in the case that Wp — 10/fci?, 
and Z = (Fig. In this case, there are several prop- 
agating modes in FMl and FM2. The magnitude of MR 
decreases with increasing L for each value of the exchange 
field. This behavior of MR is understood by considering 
the L dependence of the current in the parallel and an- 
tiparallel alignments. As shown in Fig. 0] in the case that 
ho = 0.6fiF, the finite current in the parallel alignment 



FIG. 5: MR as a function of ho for L = W/kp. 




FIG. 6: The current as a function of ho for L — lO/kp. The 
solid and dashed lines are for the current in the antiparallel 
and parallel alignments, respectively. 



flows because the ordinary Andreev reflection occurs, and 
is almost independent of L. On the other hand, the cur- 
rent in the antiparallel alignment decreases with increas- 
ing L since the contribution of the crossed Andreev re- 
flection process to the current decreases with increasing 
L, and therefore the magnitude of MR decreases with in- 
creasing L. In this case, the oscillation of the current in 
the antiparallel alignment is suppressed because electrons 
and holes in the several propagating modes I in Eq. H1U|) 
contribute to the current and wash out the oscillation. 
The reason why MR for ho = 0.8 (0.5) fip are almost 
equal to MR for ho — 0.7 (0.4) fip is as follows. In Fig. |S1 
the ho dependence of MR is plotted. We find three drops 
in MR at ho ~ 0.12 fip, 0.62 np, and 0.92 /^f- MR for 
ho = 0.8 (0.5) /ii? and ho — 0.7 (0.4) are in the same 
plateau. This plateau structure is understood by consid- 
ering the denominator /ap and the numerator Ip — Jap 
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FIG. 7: The difference between the currents in the paraUel 
and antiparallel ahgnments as a function of ho- 



FIG. 9: MR as a function of the height of the interfacial 
barriers Z for ho/fiF = 0.6. The solid, dashed, and dotted 
lines represent the case of kpL = 10, 15, and 20, respectively. 
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FIG. 8: MR as a function of L for various values of the inter- 
facial barrier parameter Z and ho/^p ~ 0.6. 



in Eq. (|29() separately. As shown in Fig. /ap is mainly 
given by the ordinary Andreev reflection, and decreases 
with increasing ho because the number of the channels for 
the minority spin decreases by one when passing across 
ho ~ 0.12 /ii?, 0.62 /ii?, and 0.92 fip. Especially, in the 
range of ho/nF — 0.92 ~ 1, there is no open channel 
for minority spin and the ordinary Andreev reflection is 
completely prohibited. Therefore, we flnd MR = — 1 (see 
Fig. ISJ. Fig. [7| shows the ho dependence of /p — /ap, 
which is mainly due to the crossed Andreev reflection. 
The magnitude of /p — /ap is much smaller than that 
of /ap , and therefore MR shows the plateau structure as 
shown in Fig. [SJ and MR for ho = 0.8 (0.5) fj^p are almost 
equal to MR for ho = 0.7 (0.4) hf- 

Finally, we investigate the effect of the interfacial bar- 
riers on the transport in this system. Figure |H1 shows the 



L dependence of MR for ho — 0.6/ii? and several values 
of interfacial barrier parameter Z . As seen in Fig. [SJ MR 
approaches zero with increasing L and shows strong de- 
pendence on the height of the interfacial barrier Z . The 
decrease of MR with increasing L is explained by the 
same way as in the case of no interfacial barriers (Fig. 
OJ. To investigate the Z dependence of MR in detail, we 
calculate the Z dependence of MR for kpL = 10, 15, and 
20 as shown in Fig. O The magnitude of MR decreases 
with increasing Z in the range of Z < 0.5 and is almost 
constant for L in the range of Z > 0.5. This dependence 
is understood as follows. MR consists of the denomina- 
tor /ap and the numerator /p — /ap , which mainly come 
from the process of the ordinary Andreev reflection and 
the crossed Andreev reflection, respectively. The crossed 
Andreev reflection is more sensitive to the scattering at 
the interfacial barriers than the ordinary Andreev reflec- 
tion, and therefore the value of /p — /ap decreases more 
rapidly than that of /ap in the range of Z < 0.5, and 
therefore the magnitude of MR decreases with increasing 
Z for kpL = 10, 15, and 20 as shown in Fig. |51 

Although the impurity scattering in SC and the prox- 
imity effect are neglected in our theory, these assump- 
tions are justifled as follows. First, as shown in the 
present calculations, the crossed Andreev reflection pro- 
cess occurs on the scale which is less than several nanome- 
ters for kp ^ lA~^m^ This scale is much smaller than 
the mean free path of SC^ and therefore the effect of 
the impurity scattering in SC on the crossed Andreev re- 
flection is neglected. Second, in the present paper, we 
consider the case that the area of the contacts of FMl 
and FM2 with SC are several nanometers and thus the 
proximity effect can be neglectedi^SiiLiSi^ 
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IV. CONCLUSION 



and 



We present a theory of the crossed Andreev reflection 
in structures consisting of a superconductor with two fer- 
romagnetic leads. By extending the BTK theory to this 
system, we calculate the current and the magnetoresis- 
tance due to the crossed Andreev reflection. It is shown 
that the dependence of the crossed Andreev reflection 
on the distance between two ferromagnetic leads, L, is 
given by the interference between the wave functions in 
ferromagnetic leads. The probability of the crossed An- 
dreev reflection follows {kpL)~^, where kp is the Fermi 
wave number, and therefore the magnetoresistance due 
to the crossed Andreev reflection strongly decreases with 
increasing kpL except for the case of half metallic fer- 
romagnets. It is also presented that the dependences 
of the magnetoresistance on the exchange field show the 
plateau structure and the magnitude of the magnetore- 
sistance rapidly decreases with increasing the height of 
the interfacial barriers. These dependences are explained 
by considering the relation between the probabilities of 
the ordinary Andreev reflection and the crossed Andreev 
reflection. 
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First, by multiplying the both sides of Eq. (|A1|I by 
^sc,miy) and integrating them with respect to y, we ob- 
tain 
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APPENDIX A: BOUNDARY CONDITIONS 

The coefficients aa,in, &<T,in, Ccr,in, d^,im ctcr,in, and Pa,ln 
in Eqs. Ijl^l) . (|14|l . and (|15|l are determined from the 
boundary conditions and as follows iSSiSSi^Si^ 
Substituting the wave functions (|13|l . (|14l) . and (|15|l for 
the boundary conditions (|16|l and H17|) . we obtain 



+E 

CO 

+E 

1=1 

oo 

= E 



*FMl,n(2/) 



1 n 



, \ 1 



*FMij(y)| 01 (y) 
*FM2,i(y) 02 (y) 
'^scAy) (y), 

(Al) 



where A2(2), ;ni(-^) is the overlap integral between the 
wave functions in FM1(FM2) and SC, and is given by 

Ai,z™(L) = / $FMi,i(2/)«'sc,m(y)d2/, (A4) 

J(L-Wf)/2 
A-L+Wf)/^ 

A2,lrn{L) = / ^FM2,liy)^SC,m{y)dy- (A5) 

J {-L-Wf)/2 



Second, by multiplying the both sides of Eq. ljA2p by 
^FMi.miy) and integrating them with respect to y, we 
obtain 
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where 6mn is a Kronecker delta defined as 



1 (m — n) 
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(A7) 



Third, by multiplying the both sides of Eq. llA2|) by 

^FM2.m{y) and integrating them with respect to y, we 
obtain 



where = p^i, = <]fi, and is the interfer- 
ence term between the wave functions of FMl and FM2 
through SC, which strongly depends on L, given by 



F±(L) = ^fc±Ai,i™(L)A2,i™(L) 



.2mH\ / 



E 

1=1 



f , 2mH\ f 



- I va 

Uo 



A. 



2,ml- 



(A8) 



By solving the Eqs. I|A3|I . (|X6)l . and (|X8)l . the co- 
efficients acr,;„, focr,i„, Ccr^in, d^,in, and ^cr,in are 
determined. In the numerical calculation, we truncate 
the number of the channels in the ferromagnetic leads 
(FMl and FM2) and SC by the cutoff constants Mp and 
Ms, respectively. The Mp and Ms are taken to be large 
enough to make the calculation results converge. 

Especially, in the half metallic {ha/iip = 1) FMl and 
FM2 with width Wf — 4/fcF, there is only one propagat- 
ing mode Z = 1 in Eq. 11U|) . In this case, we can neglect 
the evanescent mode / > 2 and take the cutoff constant 
in FMl and FM2, to be Mp = 1. In the case of no in- 
terfacial barriers {Z = 0), at low energy region {E ^ 0), 
the coefficient of the crossed Andreev reflection part in 
the wave function of FM2, c^^n, is written as 



Ci r+ + C2 r- 



(A9) 



where 



kpy/2kpWp 



ikpLy''/\xpi-L/20 



X exp ( ± i ( kpL — 



(A13) 



and 



^U2) =E^™A?(2)^i™. (A14) 
Substituting Eq. (|A9|I for Eq. lO, we obtain 



R]''II = k^Wp 



8IC3P 

X {kpLy^cxp{-L/^) 
4IC3P 

X [1 - sin {2kFL + 0)] {kpL)~^ exp (-L/C), (A15) 



where </> is a phase given by 



C2 = tp+{p- -n+){q+ + n+), 
C3 = [p+p- + n+n- -ip+ -p 



(AlO) 
(All) 

-ny/2] 
n-)/2] , 

(A12) 



and we use the relation ICiP ~ IC2P ~ ICiCn 



(A16) 
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